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ABSTRACT 

We present an analytical formula for the position and shape of the spots on the 
surface of accreting magnetized stars in cases where a star has a dipole magnetic 
field tilted at a small misalignment angle 8 < 30° about the rotational axis, and the 
magnetosphere is 2.5-5 times the radius of the star. We observed that the azimuthal 
position of the spots varies significantly when the position of the inner disc varies. In 
contrast, the polar position of the spots varies only slightly because of the compres- 
sion of the magnetosphere. The azimuthal width of the spots strongly varies with 8: 
spots have the shape of an arc at larger misalignment angles, and resemble a ring 
at very small misalignment angles. The polar width of the spots varies only slightly 
with changes in parameters. The motion of the spots in the azimuthal direction can 
provide phase-shifts in accreting millisecond pulsars, and the "drift" of the period in 
Classical T Tauri stars. The position and shape of the spots are determined by three pa- 
rameters: misalignment angle 8; normalized corotation radius, r c /R* and normalized 
magnetospheric radius, r m /R*, where R* is the stellar radius. 

We also use our data to check the formula for the Alfven radius, where the main 
dependencies are r m - (|li 2 /]VL) 2/7 , where ill is the magnetic moment of the star, and M 
is the accretion rate. We found that the dependence is more gradual, r m - (|li 2 /]VL) 1/5 , 
which can be explained by the compression of the magnetosphere by the disc matter 
and by the non-dipole shape of the magnetic field lines of the external magnetosphere. 
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1 INTRODUCTION 

Magnetospheric accretion occurs in a variety of astrophysical 
systems. The accreting matter is stopped by the stellar mag- 
netic field roughly at a distance from the star where matter 
and magnetic stresses become equal. Beyond that point, mat- 
ter flows around the magnetosphere in a funnel flow and falls 
near the magnetic poles of the star, forming hot spots (e.g., 
|Ghosh&Lamb|1979| ). 

The light-curves observed from accreting magnetized 
stars are often associated with the hot spots on their surfaces. 
In many applications it is important to know the exact loca- 
tion and shape of the spots. For example, in millisecond pul- 
sars, pulse profiles are significantly affected by the hot spot 
location and shape (e.g., Poutanen & Gierlinski 2003; Ibragi- 



|mov & Poutanen|2009| |Leahy et al.|2009p . The position and 
shape of the spots is also important for understanding the 
light-curves of Classical T Tauri stars (CTTSs, e.g., |Bouvier| 



et al.|2007|) a nd in magnetized accreting white dwarfs (e.g., 
Warne r 1 19951 ). 
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Three-dimensional (3D) magneto-hydrodynamic (MHD) 
simulations of magnetospheric accretion show that matter ac- 
cretes to a star in two funnel streams, which form two spots 
on its surface dRomanova et al.||20"03] |2004l |Kulkarni & Ro] 
manova 2005). Simulations show that the spots have the 
shape of a ring at small misalignment angles, the shape of an 
arc at 8 « 30° and the shape of a bar at very large 8. In all 
the cases, the energy flux is largest in the center of the spot 
and gradually decreases outward. The position of the spots 
does not coincide with the magnetic pole. 

The simulations also show that the funnel stream can be 
dragged by the disc so that the spot forms at higher longi- 
tudes on the star's surface. Alternatively, it can trail the mag- 
netosphere, causing the spot to form at lower longitudes ( |Ro-| 
man ova et al.|2004] ) . Therefore, the azimuthal position varies 
with accretion rate. In the case of a small dipole misalign- 
ment angle, 8 < 5°, the funnel stream may pass a whole 
cycle about the magnetic pole, so that a spot can move faster 
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or slower than the star (e.g., Kulkarni & Romanova|2008[|Ba-| 
Ichetti et al.|2010| ). 

The motion of the spots along the stellar surface can 
produce observable effects such as the phase shifts in light- 
curves, timing noise and intermittency of accreting millisec- 
ond pulsars (e.g.,|Lamb et al.|2008||2009||Patruno et al.|2009l 
|Poutanen et al.||2009| ). The motion of the spots can possibly 
explain the drifting periods observed in many CTTSs (e.g., 
Rucinski et al. 2008 ). Usually a simple model is used for the 
spots, such as a circular spot centered at the magnetic pole, 
with a constant or gaussian distribution of emitted flux. 

More recent numerical simulations show that accreting 
magnetized stars may also be in an unstable regime of ac- 
cretion, where matter accretes due to the magnetic Rayleigh- 
Taylor instability (Kulkarni & Romanova 2008; Romanova 
et al. 2008]). In this regime, multiple chaotic 'tongues' pen- 
etrate through the magnetosphere and form spots of chaotic 
shapes and positions. This regime is favorable when the star 
rotates slower than the inner disc, though additional factors 
are important for the onset of instability (e.g., | Spruit et al.| 
1995 ). In this paper, we only consider the set of parameters 
at which accretion is stable and the spots have regular shapes. 

The goal of this work is to derive a convenient formula 
for the shape and position of the spots, and the dependen- 
cies of the parameters of the spots on the parameters of the 
star and the disc. We concentrate on the cases of relatively 
small misalignment angles, 6 < 30°, where the spots have 
the shape of an arc or a ring. 

We also use a set of our 3D MHD simulations to test the 
standard formula for the Alfven radius, which is routinely 
used for calculation of the truncation (magnetospheric) ra- 
dius, r m . We found that the dependencies of r m on the mag- 
netic moment of the star, \i, and the accretion rate, M, are 
different compared with the standard formula. 

In Sec. 2 we describe our numerical model. In Sec. 3 we 
provide an analytical formula for the spots. In Sec. 4 we in- 
vestigate the dependence of the parameters of the spots on 
the parameters of the star and the disc. In Sec. 5 we provide 
a brief practical guide on how to find the position of the spots 
from observation. In Sec. 6 we use our simulations to test the 
formula for the Alfven radius. In Sec. 7 we outline the main 
conclusions from this work. Appendix A describes our refer- 
ence values. Appendix B shows parameters of the spots for 
all simulation runs. In Appendix C we derive the formula for 
the magnetospheric radius and compare it with the standard 
formula for the Alfven radius. 



2 NUMERICAL MODEL 

To obtain the shapes of the spots, we perform global three- 
dimensional (3D) magnetohydrodynamic (MHD) simulations 
of accretion from the disc onto a magnetized rotating star. 

The model we use is the same as in our earlier 3D 
MHD simulatio ns (e.g., |Romanova et al.|2004[|Kulkarni & Ro-| 
manova 2008 ). The star has a dipole magnetic field, the axis 
of which makes angle 6 with the star's rotational axis. The 
rotation axes of the star and the accretion disk are aligned. 
A magnetized star is surrounded by an accretion disc which 
is cold and dense, and by a hot rarefied corona which is 
100 times less dense and 100 times hotter in the fiducial 
point. The disk and the corona are chosen to be initially in 



a quasi-equilibrium state, where the gravitational, centrifu- 
gal and pressure gradient forces are in balance (Romanova 
|et al. ||2002] ). General relativistic effects, which are impor- 
tant in neutron stars, are modelled using the Paczyhski-Wiita 
potential, <D(r) = GM,/(r-r g ) ( [Paczyhski & Wiita|[l980] ), 
where M* is the mass of the star and r g = 2GM*/c z is 
Schwar zschild radiuj*] Viscosity is modelled using the oc- 
model ( [Shakura & Sunyaev|1973[|Novikov & Thorne|1973| ). 
It is incorporated only in the disc and controls the accretion 
rate through the disk. We take a small parameter oc = 0.02 in 
all simulation runs. 

To model accretion, the MHD equations are solved nu- 
merically in three dimensions using a Godunov-type numeri- 
cal code, written in a "cubed-sphere" coordinate system which 
rotates with the star ( Koldo ba"et al.|2002]). The nume rical ap- 
proach is similar to that described in Powell et al.| ([1999 ), 
where the eight-wave Roe-type approximate Riemann solver 
is used to calculate flux densities between the cells. The grid 
resolution is identical to that in Kulkarni & Romanova ( 2008 ) 
and equal toN r x N 2 = 72 x 31 2 in each of the six blocks 
of the cubed sphere. Here, N r is the number of grids in the 
radial direction, and N is the number of grids in the angular 
directions in each of the six sides of the cube. 

The boundary conditions at the star's surface amount to 
the assumption that the infalling matter passes through the 
surface of the star, so that the dynamics of this matter after it 
falls onto the star is ignored. At the external boundary, mat- 
ter inflow is permitted at the disc's part of the boundary, and 
outflow at the corona's part of the boundary. 

The simulations are done using dimensionless variables. 
The dimensionless value of every physical quantity q is de- 
fined as q = q/qo, where q is the reference value for q. Ap- 
pendix^ shows how the reference values are determined, and 
lists the reference values for three classes of central objects: 
classical T Tauri stars, white dwarfs and neutron stars. Subse- 
quently, we drop the tildes above the dimensionless variables 
and show dimensionless values everywhere unless otherwise 
stated. For clarity, we keep the tildes above the dimensionless 
magnetic moment of the star, and the matter and angular 

momentum fluxes M and L, respectively. 

Fig. [T] (left panel) shows an example of funnel streams 
obtained in a typical 3D MHD simulation run. The right panel 
shows the energy flux distribution at the surface of the star 
(the hot spot) . One can see that the spot has the shape of an 
arc, and the center of the spot (the highest energy flux) is 
off-set from the magnetic pole. The center of the spot is also 
off-set from the u. — Q. plane, because the inner disc rotates 
faster than the magnetosphere, and both the funnel stream 
and the spot are shifted. In this paper, we provide an analyti- 
cal formula for the shapes and positions of such spots. 



1 Since the simulations include general relativistic effects (approxi- 
mated using the Paczyriski-Wiita potential), strictly speaking they are 
only applicable to neutron stars. However, previous studies ( Kulkarni 
& Romanova 2005) show that use of the Paczyriski-Wiita potential 
does not affect the shape of the spots significantly. So the results in 
this work can be applied to other types of accreting systems as well. 
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Figure 1. Left panel: A 3D view of the funnel flow from the disc to a magnetized star, where the dipole moment \x is tilted by = 20° about the 
rotational axis. One of the density levels is shown in green; sample field lines are shown in red. Right panel: the energy flux distribution on the 
surface of the star. Circles show the position of the magnetic and rotational (O) axes, respectively. Other parameters are ft = 2 and r c = 2. 
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Figure 2. Left panel: Energy flux distribution on the surface of the star obtained in one of the simulation runs, with parameters = 10°, jl = 1.5, 
r c = 1.8, at time t = 8. Middle panel: Normalized flux distribution in the fitted spot, calculated with Eq.[T] The white line shows the boundary 
where the flux is e— times smaller than in the center of the spot. The position of the spot in polar, O , and azimuthal, cj) > directions is shown. 
Right panel: The same as the middle panel, except that the width of the spot is denned in polar, 6i, and azimuthal, 4>i, directions. 



3 ANALYTICAL FORMULA FOR THE SPOTS 

For relatively small misalignment angles, < 30°, the posi- 
tion and shape of the spots can be approximated by analytical 
formula. Namely, each antipodal spot can be well approxi- 
mated as a circular arc centered at the magnetic pole, with a 
gaussian flux distribution centered some distance away from 
the magnetic pole, as follows (see Fig. [2] two right panels) : 



F(9, 4>) = F c exp 



• 9 V , {§-(bo x 2 



9i 



(1) 



The spherical polar and azimuthal angles 9 and are mea- 
sured with respect to the magnetic axis \i. The azimuthal posi- 
tion (\) = is defined by the \i— D. plane. The spot is described 
by the following parameters: 

• The polar and azimuthal position of the spot (9 , <t> ). 

• The polar and azimuthal width of the spot (9 la cj>i), which 
are determined to be a half-width of the gaussian in the polar 
and azimuthal directions. The half- width is defined as the dis- 



tance from the center of the spot, to the radius at which the 
emitted flux is 1/e times smaller than at the spot center. 
• The flux emitted from the center of the spot, F c . 



We then apply a least-squares fit of this expression (given 
by Eq. [TJ to the spots obtained from the MHD simulations. 
Fig. [5] compares the northern hot spot from the simulation, 
where misalignment angle 6 = 10° (left panel), with the fit- 
ted spot (right two panels). The fitted spot agrees well with 
the simulation. This procedure, when applied to the spots ob- 
tained at different parameters, shows that the rms error of the 
approximation is usually a few per cent. 

Equation [T] gives the flux distribution for the northern 
hot spot. The southern hot spot is generally identical to the 
northern one for a dipole field, so that parameters describ- 
ing it are simply [n — 9 0j 7t + cj) ) and (9 X , cj)i). When using 
Eq. [T] with these parameters, care must be taken to ensure 
that \§ — (f>ol < 7t (by adding integral multiples of n to <£) for 
northern spots and I4> — (7t + (J> )| ^ n for southern spots. 
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4 DEPENDENCE OF SPOT PARAMETERS ON 
PARAMETERS OF THE STAR 

We performed a large set of 3D MHD numerical simulations 
at different misalignment angles of the dipole moment 8, pe- 
riods of the star P* and its magnetic moment \L The disc is 
the same in all simulation runs. Our simulations and analysis 
show that the location and shape of the spots are determined 
by the following parameters: 

(i) Misalignment angle 8. We varied this parameter from 
very small values up to relatively large values: 8 = 
2°, 5°, 10°, 20°, 30°. 

(ii) Magnetospheric radius r m , normalized by stellar radius 
r m /R*. The magnetospheric radius can be determined in sev- 
eral ways. For example, an equality of matter and magnetic 
stresses, |3i = (p + pv^)/(B 2 /87t) = 1, usually corresponds 
to the innermost region of the disc, where the disc density 
decreases sharply (e.g., Romano va et al.|2011| ). On the other 
hand, the condition (3 = 87ip/B 2 often corresponds to the re- 
gion in the disc where the funnel flow begins ( Bess olaz et al.| 
2008 ). Our simulations show that the condition (3 = 1 corre- 
sponds to the radius from which matter flows along a funnel 
stream to the central regions of the spots. This is why we 
chose the condition (3 = 1 to derive r m . We varied the size 
of the magnetosphere using the dimensionless magnetic mo- 
ment of the star ft in the range of ft = 0.5 — 2, and obtained 
the magnetospheric radius in the range of r m /R* « 2.4 — 4.9. 

(iii) Dimensionless period of the star P*, or the ratio r c /R*, 
where r c = [GlVU(P*/27t) 2 ] 1/3 is the corotation radius (here 
P* is the dimensional period). In our model, we use either 
dimensionless period P*, or the ratio r c /R*. We chose three 
periods for the star: P* = 1.8,2.4,2.8 which correspond to 
r c /R* = 4.2, 5.1, 5.3 (see Table [T] for the conversion between 
dimensional and dimensionless values) . 

For each simulation run, we chose an interval of time 
where the accretion rate is approximately constant. We de- 
rived an averaged spot using several moments in time from 
this interval, then approximated the observed spot (energy 
flux distribution) with Eq. [T] and derived best-fit parameters 
(0 O , cf>oL (9i, and F c . Table [bT] lists these parame ters for 
different values of 8, r m /R*, and P* (r c /RJ. Figures 



Bl 



B2 



and [B3] show each of the five parameters of the spot for all 
simulation runs (symbols) for different misalignment angles 
8 (different colors) and different periods of the star P* (dif- 
ferent symbols) . Below, we discuss each parameter and derive 
useful dependencies. 

4.1 Polar position of the spot, 9 

We derive the polar position of the spot using two methods: 
(1) simple analytical approach for a pure dipole field, and (2) 
numerical simulations. 

4.1.1 Polar position of the spot derived from the dipole model 

One of the parameters - the polar position of the spot, 9 , - 
can be derived analytically. 

Fig. [3] (left panel) shows a schematic view of a star with 
a dipole field tilted at angle 8. The dipole field line which 
crosses the inner disc radius r = r m is shown. This line also 
shows the expected position of the spot. The dipole field lines 



obey the relation r/ sin 2 9 = constant, where 9 is measured 
from the magnetic axis. For simplicity, let us assume that the 
center of the spot is in the \i — D. plane. 
Referring to Fig. [3] we have 



(2) 



8 



sin 2 9/ 



which gives us the location of the spot, 9 S , as 



e s 



: sin 




cos 8 



(3) 



In this idealized model, the derived angle 9 S shows the 
approximate position of the spot in the polar direction, 9 ~ 
9 S . The position of the spot depends on truncation radius 
r m /R* and misalignment angle 8. 

This could be a useful approach in finding the polar po- 
sition of the spots. However, numerical simulations show that 
the polar position is different from that derived analytically, 
because the magnetic field is compressed by the disc. Com- 
pression changes the shape of the field lines, affecting the po- 
sition of the spots on the star. 

4. 1 .2 Polar position of the spots derived from numerical 
simulations 

Table |bT] and Fig. [Bl] (left panel) show values of 9 for all 
simulation runs. It can be seen that the polar position varies 
within a narrow range of angles, 16° < 9 < 18°. No any sig- 
nificant correlation has been found between 9 and the other 
parameters. This result is somewhat unexpected, because ac- 
cording to the simple estimates performed in the previous sec- 
tion, larger values of 9 are expected for smaller misalignment 
angles 8 and smaller magnetospheres (smaller r m /R*), and 
vice-versa. According to Eq. [3] if we use extreme parameters 
of 8 = 2° (the smallest misalignment angle) and r m /R* = 2.4 
(the smallest magnetosphere), then we expect the largest an- 
gle to be 9 = 40.2°. For 8 = 30° and the largest magne- 
tosphere, r m /R* = 4.9, we expect the smallest angle to be 
9 = 23°. Both angles are larger than the angles obtained in 
our numerical simulations. 

To understand this phenomenon, we compared the 
shapes of the magnetospheres obtained in our simulations 
with the shape of the dipole magnetosphere. An example of 
such a comparison is shown in Fig. [3] (right panel). Within the 
inner magnetosphere, the field lines of the simulated magne- 
tosphere (yellow lines) depart only slightly from the dipole 
field line (shown in red). However, in the outer magneto- 
sphere, the field lines of the modeled magnetosphere strongly 
depart from the dipole field lines. Matter flowing along the 
compressed field line forms a spot closer to the magnetic pole 
(at smaller 9 ), compared with the theoretically-predicted po- 
sition of the spot in case of a dipole field. 

Fig. [3] also shows the position of the magnetospheric ra- 
dius r m , where the (3 = 1 line crosses the equatorial plane 
(see two black circles) . It can be seen that the field lines which 
start at r m connect the disc with the central parts of the spots. 

We also plot slices of density distribution in several cases 
with "extreme" parameters. Namely, we choose those simula- 
tion runs where parameters 8 and r m /R* are very large or 
very small. Thus, according to Eq. [3] the spots are expected 
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Figure 3. Left Panel: The sketch shows a star with a dipole magnetic field, the position of the inner disc r = r m and the polar angle 9 S of the 
spot's position on the star. The magnetic axis \jl of the dipole moment is tilted at an angle about the rotational axis O.. Right Panel: Density 
distribution (color background) and magnetic field lines (yellow lines) in a typical simulation run (0 = 30°, r c = 1.8, jl = 1.5, at t = 15). 
Red lines show the dipole field lines at t = 0. The external field lines are chosen such that both the yellow and the red lines start at the inner 
magnetospheric radius r m , which is marked as a black circle. 
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Figure 4. The plots show that the polar angle 9o (the angle between solid lines) is approximately the same for different parameters of the star. 
Top and bottom panels show cases of large (0 = 30°) and small (0 = 5°) misalignment angles. The left and right panels show cases of small 
(jl = 0.5) and large (jl = 1.8) magnetic moments of the star. Light-gray line shows The (3 = 1 line is shown in light gray, and the yellow dots 
show the position of the magnetospheric radius r m . The dipole field lines, which cross the magnetospheric radius are shown in cyan. The dashed 
black lines show the theoretically-predicted polar position of the spot, 9 S (see Eq.|3}. 



to be at different polar positions 9 S . Fig. [4] shows that in spite 
of the difference in parameters, the polar position of the spot 
9 is approximately the same and corresponds to the narrow 
range of values observed in Fig. [M] This can be explained 
by the fact that in cases of a large magnetosphere compres- 
sion is not very strong, but is sufficient enough to shift the 
spot position towards the magnetic pole and yield a value 
of 9 which is different from that predicted theoretically in 
the case of a dipole field. In cases of small magnetospheres, 
theoretically-predicted angle 9 S is larger than in cases of a 
large magnetosphere. However, compression of the magneto- 



sphere is stronger, and the spot has a stronger shift towards 
the magnetic pole, compared with cases of a large magneto- 
sphere. This effect leads to similar angles 9 in our simulation 
set. It is clear that at even smaller (than in our set) magneto- 
spheres, the theoretical dependence will not work, due to an 
even stronger compression of the magnetosphere. However, 
the formula may be applicable at very large r m /R*, where 
compression is expected to be weaker. 
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Figure 5. Left Panel: Dependence of 4>o on the magnetospheric radius r m /R* for different periods of the star P*. The case of 8 = 2° is removed 
from the set because the spot changes its position, and its averaged preferred position has a large error. Right panel: Dependence of the azimuthal 
position of the spot cf>o on the ratio of magnetopsheric to corotation radius r m /r c , taken for all simulation runs. 



4.2 Azimuthal position of the spot, cj) 

Figure |bT] (right panel) shows the azimuthal position of the 
spots for all simulation runs. The figure shows a wide range 
of angles cj) . We noticed that some correlation can be seen 
between groups of runs with the same period of the star P* 
(which are marked with the same symbol in the figure). To 
better demonstrate this apparent correlation, we divided all 
data into three groups corresponding to three periods of the 
star, and obtained much clearer dependencies of cj) on r m /R*. 
The list-square approximation shows linear correlations: 



cj)o = 84° - 17° 



for P* 



1.8 



4.3 



(4) 



* = 206° - 39° 



for P* = 2.4 



* = 300° 



57° -p, for P* 



2.8 



R* 



5.1 (5) 



5.7 (6) 



It can be seen that angle <$> varies more rapidly in the 
cases of a slower- rotating star (P* = 2.8), and more slowly 
in cases of a more rapidly- rotating star (P* = 1.8). 

These dependencies can be interpreted as an azimuthal 
phase-shift that varies with accretion rate: smaller values 
of r m /R* correspond to higher accretion rates and stronger 
shifts cj) . 

The azimuthal position of the spot depends on the rel- 
ative velocity of the inner disc with respect to the magneto- 
sphere. Hence, the dependence of <fy on the ratio r m /r c is 
expected. In this case, we were able to obtain a unique corre- 
lation for all simulation runs. Fig. [5] (right panel) shows this 
correlation. The least-square fit shows a linear dependence: 



4>o = 160° - 146° 



(7) 



Therefore, if a value of r m /r c is known, then the approx- 
imate azimuthal position of the spot is known as well. 

Fig. [^demonstrates how the spot position cj) varies with 
^m/R* (which can be interpreted as variation of the accretion 
rate): angle (j) increases when r m /R* decreases. 

We should note that in the case of very small misalign- 
ment angle, 8 = 2°, the spot changes its position rapidly, 



making entire 2?t cycle around the magnetic pole. Hence, the 
azimuthal position of the spot has a large error. This is why 
we excluded the case of 8 = 2° from the left panel of Fig. [5] 



4.3 Polar width of the spot, 9i 

Figure [B2] (left panel) shows that the polar width of the spots 
varies within a relatively narrow interval, 7° < 0i < 10°, and 
is Gi = 7° — 8° for most models. In a few cases, it is larger: 
0i = 9°, and 9 = 10° in one simulation run. There is a weak 
correlation between 0! and r m /R*, with 0! being smaller in 
the cases of larger r m /R*. 



4.4 Azimuthal width of the spot, <$>i 

Figure |B2] (right panel) shows the azimuthal width of the 
spots fyi for all simulation runs. This parameter is also scat- 
tered significantly, like <$> . However, some ordered varia- 
tion can be seen if we consider the subsets of runs with the 
same misalignment angle 8 (marked by different colors) . The 
points with a large 8 are located at the bottom part of the 
plot, where fyi « 100° and the spots have the shape of an 
arc. The points with smallest 8 are located towards larger 
<$>i, and have the shape of a ring. 

Taking into account this (expected) dependence of c()i on 
8, we separate all the runs into subsets with a different 8. We 
notice that the correlations become stronger when we plot c()i 
as a function of r m /r c . Fig. [7] (left panel) shows subsets of 
data for each 8 and a linear approximation for each subset. 
At large misalignment angles, 8 = 20° — 30°, the azimuthal 
width does not vary much with r m /r c , and fyi « 90° — 100°. 
However, c()i strongly increases with r m /r c for small 8. The 
lines which approximate sets of data have different slopes, 
and thus, the dependencies on 8 are different. To derive these 
dependencies, we choose several values of r m /r c and take 
the values of $i that correspond to lines at different 8. We 
plot these values in Fig. [7] (right panel) and also obtain the 
following dependencies: 
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Figure 6. The plot demonstrates different azimuthal positions (£o of the spots at different magnetospheric radii r m /R*. The spot shifts in the 
azimuthal direction when r m /R* decreases. The colored background shows the energy flux, which varies from the minimum value F m i n = 0.5 
(dark cyan color in all panels) to the maximum value F max (red color), which is 5.8, 4.8, 3.2 and 1.9, from left to right. 



*1 


= i50°e- 


-0.03 


for 


r m /r c 


= 0.6, 


(8) 


*1 


= 400°e- 


-0.37 
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for 


r m /r c 


= 0.8, 


(9) 


*1 


= 730°B- 


-0.60 

3 


for 


r m /r c 


= 1.0 . 


(10) 



The Figure and Eqs. 8-10 show that for 8 > 15°, the 
width of the spots is relatively small for all values of r m /r c . 
However, for 8 < 15°, the shape varies from an arc for 
smaller r m /r c to a ring for larger r m /r c . For a star with the 
same corotation radius r c , larger/smaller values of r m /r c cor- 
respond to a smaller/larger accretion rate. 

Fig. [8] demonstrates the shape of the spots in the case 
of the same magnetic moment jl and period P* of the star 
(approximately the same value of r m /r c ), but a different mis- 
alignment angle of the dipole 8. As we can see, the shapes of 
the spots vary from a ring (at small 8) to an arc (at large 8) . 

4.5 Energy flux in the center of the spots, F c 

From our simulations and fits of the observed spots with Eq.[T] 
we obtain the maximum energy flux in the center of the spot 
F c (in dimensionless units). The dimensional flux F c , which 
is used in Eq. [I] can be derived from Eq. |A2| It can be seen 



that the main dimensionless variable in Eq. is the ratio 



Bl 



F c = F C /]I 2 . This is why we put this variable into the Table 
and use it in Fig.|B3| 

Fig. B3 shows that the values of F c are mainly confined 
to the interval 0.5 < F c < 4. It is larger, F c « 6.5 - 8.5, 
in a few points corresponding to a smaller magnetosphere, 
r m /R + « 2.8 — 3.2. The fluxes are about two times smaller in 
cases of a small 8 (see red line in Fig. |B3]) . 



5 MAGNETOSPHERIC RADIUS 

Magnetospheric radius r m (which corresponds to the bound- 
ary between the magnetosphere and the disc) is often derived 
from the condition that it is proportional to the Alfven radius 
r° A , derived for spherical accretion under the assumption that 
the magnetic pressure in the magnetosphere equals to the ram 
pressure of free-falling matter (e.g., |Lamb et al.|1973| ): 



V7 



(ID 



V 2GM*M. 2 , 

where k A is the dimensionless coefficient of the order of unity. 
This formula has been tested a few times. In these tests, the 



inner disc radius has been derived in two different ways. First, 
it has been derived from the simulations using conditions 
3 = 1 or |3i = 1. Second, it has been derived using formula 
11 for r° A , where values such as \i and M were taken from the 



simulations. |Long et aL] ( |2005| ) performed such comparisons 
in three models with slightly different parameters u. and the 
condition pi = 1, and obtained an average approximate value 
of k A w 0.5. [Bessolaz et al.| ( |2008| ) used a variety of defini- 
tions for the inner disc radius, such as (3 = 1 and |3i = 1, and 
obtained k A ~ 0.4 in the case of the (3 = 1 condition (see also 
Zanni & Ferreira 2013 ). However, none of this work studied 
the dependence of r m on \i or M. 

We have a sufficiently large set of simulations, so these 
dependencies can be tested here. As the first step, we take the 
formula for the Alfven radius (Eq. [TT] ), substitute in our di- 
mensional units (see Sec.|A|), and derive the magnetospheric 
radius in a dimensionless form: 



k A 
0.35 



2/7 



(12) 



We can test this formula in our simulations, where the di- 
mensionless magnetic moment jl is the main initial parameter 
in each simulation run, and the dimensionless accretion rate 

M is the output from each simulation run. We also derived 
the ratio r m /R* from the condition that |3 = 1 in each simula- 
tion run. Fig. [9] shows the resulting dependence of r m /R* on 

p 2 /M. The dependence is well-approximated by the power 
law in the following form: 



i 2.2 



~2\ 1/5 

M 



(13) 



We can see that the power obtained from the simulations, 
a « 0.2 = 1/5 is smaller than the power given by the for- 
mula for the Alfven radius, a = 2/7 ~ 0.29. We suggest that 
the compression of the magnetosphere is probably responsi- 
ble for different dependence of r m on (I 2 /M.. Namely, when 
the magnetosphere is compressed, the magnetospheric radius 
varies slower with magnetic moment and accretion rate. 

From the comparison of equations [T3] and [T2] we can de- 
rive coefficient k A , which is not a constant anymore: 



(14) 
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0.8 /r 0.9 
m c 

Figure 7. Left Panel: Dependence of the azimuthal width of the spot cf>i on the ratio of magnetic to corotation radius r m /r c for different 
misalignment angles 0. Right Panel: Dependence of (£1 on the misalignment angle of the dipole moment for different r m /r c . 






0=10° 




©=20° 




©=30° 




Figure 8. Shapes of the spots in cases of different 0. The colored background shows the energy flux distribution. The black lines indicate the 
boundary where the energy flux is e times smaller than the maximum flux. Other parameters are: ft = 1, r c = 1.8, and are the same in all cases. 
The energy levels vary from the minimum value of F m i n = 0.5 (light-cyan color), which is the same in all cases, to the maximum value of P max 
(red color), which is 2.1 in the three left panels, and is 3.0 and 3.8 in the 4th and 5th panels, respectively. 



The dependence of k A on \± 2 /M is weak, so that for the 

range of2<|I 2 /M<50 used in our set, we obtain k A in the 
range of 0.55 < k A < 0.72. These values of parameter k A are 
in agreement with the values derived earlier in axisymmetric 
simulations (e.g., |Long et al.|2005] ). 

Next, we use Eq.[T3]and Eq.|XT]to derive the ratio r m /R* 



through dimensional parameters: 



R* 



1.06 



M 2 G1VUR7 



1/10 



(15) 



This formula can be used if the magnetic moment of the star 
and the other parameters are known. 

We should note that the dependence derived in Eq.[T3]is 
relevant for the magnetospheric radii 2.5 < r m /R* < 5 inves- 
tigated in this work. A separate analysis should be done for 
larger and smaller radii. We anticipate that at larger r m /R* 
compression will be less important and the power in Eq. \Ts\ 
will be larger than 0.2, while at smaller r m /R* compression 
will be more important and the power will be smaller. Spe- 
cial simulation runs should be performed to investigate this 
dependence at larger and smaller r m /R*. 



5 

4.5 



2.5 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



(r/R.) = 2.2 (£ 2 /M) 1/5 




30 40 50 



Figure 9. Normalized magnetospheric radius r m /R* as a function of 

the ratio ft 2 / M for all simulation runs (circles) . The line shows the 
power fit approximation. 

6 SOME PRACTICAL CONSIDERATIONS 

The above analysis shows that the shapes and positions of the 
spots are determined by three parameters: (1) the misalign- 
ment angle of the dipole moment 0, (2) the normalized coro- 
tation radius r c /R* and (3) the normalized magnetospheric 
radius r m /R*. Observations of magnetized stars are usually 
complex, and most of the parameters of the star and the inner 
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P* (dim) 


r c 


r c /R* 


P* (CTTSs) 


P* (White dwarfs) 


P* (Neutron stars) 


1.8 


1.5 


4.3 


3.2 days 


52 s 


4.0 ms 


2.4 


1.8 


5.1 


4.3 days 


70 s 


5.3 ms 


2.8 


2.0 


5.7 


5.0 days 


81 s 


6.2 ms 



Table 1. From left to right: dimensionless period of the star P*, dimensionless corotation radius r c ; corotation radius normalized to stellar radius, 
r c /R*; dimensional period of the star P*, for different types of stars. 



disc are not known, or can be estimated only approximately. 
Below, we discuss possible ways to derive these parameters. 

• It is often the case that the period of the star P* is one 
of the best-known parameters. If mass and radius R* of 
the star are known, then the dimensionless period P* and the 
ratio r c /R* can be derived. 

• Many evolved magnetized stars are expected to be in the 
rotational equilibrium state, where a star's spin varies around 
the value corresponding to rotational equilibrium. In these 
cases, the magnetospheric radius can be derived from condi- 
tion r m = kr c (where k « 1) If a star is not in rotational 
equilibrium, then this condition cannot be used. 

In another approach, if the magnetic moment of the star 
and the accretion rate are known, then r m can be calculated 
from the formula for the Alfven radius derived for spherical 
accretion (see Eq. [TT] ) . In this paper, we perform an indepen- 
dent derivation of a similar formula and find that it differs 
from the standard formula for the Alfven radius (see Sec. [5}. 

• One of the most important parameters, fa, which stands 
for the azimuthal shift of the spot, depends only on the ra- 
tio r m /r c (see Eq. [5]). If a star is in the rotational equilib- 
rium state, then the value r m /r c « can be taken from nu- 
merical simulations. Axisymmetric simulations performed by 
Long et al.|([2005| ) show that r m /r c « 0.7 - 0.8, while |Zann"I| 
8c Ferreira ( 2013| ) found r m /r c « 0.6. This ratio is smaller, 
when a larger amount of angular momentum flows into the 
corona along inflated or partially-inflated field lines. This ra- 
tio can also vary in time. Therefore, the above numbers are 
approximate and can be used only as estimates. 

• The misalignment angle of the dipole's magnetic moment 
0, is known approximately in some cases. For example, in a 
few CTTSs, the magnetic field distribution on the surface of 
the star has been derived from polarimetric observations and 
then approximated with the set of tilted magnetic moments 
of a different order. It was found that the misalignment angle 
of the dipole component is usually small, 6 « 10° — 20° (e.g., 
(Donati et al. 2008)). In some cases, the misalignment an- 
gle can be obtained from the shapes of observed light-curves 
and comparisons with those obtained in 3D MHD numerical 
simulations ( [Romanova et a l. 2004). For example, for some 
probable inclination angle of i « 45°, a nearly- symmetric 
light-curves with one peak per period yield a relatively small 
misalignment angles, 6 < 30°, while light-curves with two 
peaks per period yield a larger 6. 



2 The condition k = 1 corresponds to the case where the magneto- 
sphere interacts only with the innermost radius of the disc at r = r m . 
In reality, stellar angular momentum may partially flow to the corona 
along inflated field lines, and partially to the disc along the field lines 
connecting a star with the disk at r > r m . Therefore, it is expected 
that k < 1 in rotational equilibrium state. 



7 CONCLUSIONS AND DISCUSSION 

In this work we derived a useful analytical formula (Eq. 
[TJ that provides an approximation for arc-shaped and ring- 
shaped spots which form in accreting magnetized stars. Mul- 
tiple 3D MHD simulation runs were performed to find the 
positions of the spots (9 , <$>o), their shapes (9i, <$>{) and their 
dependencies on the period of the star P*, the misalignment 
angle of the dipole 6, and the ratio of the magnetospheric 
radius to the radius of the star, r m /R*. The main conclusions 
are as follows: 

1. The polar positions of the spots 9 vary within the nar- 
row range of 16° < 9 < 18°. These angles are different from 
the polar angles obtained with the simple analytical formula 
(see Eq. [3]), where 9 has larger values and varies consider- 
ably with r m /R* and the dipole misalignment angle 6. We 
noticed that in all our numerical models the external closed 
magnetosphere is compressed and departs significantly from 
the dipole shape, which leads to the misplacement of the spot. 

2. The azimuthal positions of the spots fa strongly vary 
with the magnetospheric radius r m /R* (accretion rate). They 
increase, when r m decreases. For the same interval of r m /R*, 
fa varies more rapidly for stars with larger periods. 

3. An important correlation has been found for the de- 
pendence of fa on the ratio r m /r c , which is valid for all sim- 
ulation runs. 

4. The polar widths of the spots 9i vary within the nar- 
row interval of 7° < 9i < 9°, and the majority of the spots 
have a width of 7° -8°. 

5. The azimuthal widths of the spots $1 strongly depend 
on the dipole misalignment angle 8. They vary from fyi = 
90° - 100° for 8 = 30° (arc-shaped spots) to fa « 600° for 
Q = 2° (ring-shaped spots) . 

6. We also used our data to check the formula for 
the Alfven radius, where the main dependence is r m - 
{\i 2 /M) 2/? . We found that the dependence is more gradual, 
r m - (\± 2 /M) 1/s , which can be explained by the compres- 
sion of the magnetosphere by the disc matter and by the non- 
dipole shape of the magnetic field lines of the external mag- 
netosphere. 

Once these parameters are known, the analytical formula 
given by Eq. [I] can be used to model the energy flux distri- 
bution in the spot. The light-curves can be calculated using 
a separate numerical program. In this program, the radiation 
model, the anisotropy of the radiation, and the inclination an- 
gle of the system can be taken into account. For neutron stars, 
this code should take into account relativistic light-bending, 
Doppler-shift and other relativistic effects (see, e.g., Kulkarni 
|&Romanova|2005| ). 

Variation of the azimuthal position of the spot fa with 
accretion rate is particularly important for observations. For 
example, in millisecond pulsars, it may lead to phase-shifts 
in light-curves ([Patruno et al.||2009| ) and timing noise (e.g., 
|Poutanen et al.||2009| ). In CTTSs, it can lead to the phe- 
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nomenon of a drifting period, where the period derived at 
different times of observation, is different, and therefore, it 
varies with time (e.g., |Rucinski et al.|2008| ). 

Here, we present results for a relatively narrow inter- 
val of the ratio 0.6 < r m /r c < 1.1, which corresponds to 
ordered accretion through two funnel streams. The smallest 
values correspond to a boundary between stable and unstable 
regimes of accretion^] where at r m /r c < 0.6 accretion pro- 
ceeds in the equatorial plane through the magnetic Rayleigh- 
Taylor instability. The largest value, r m /r c « 1.1, corresponds 
to the propeller regime, in which the funnel accretion is for- 
bidden by the centrifugal barrier (e.g., |Illarionov & Sun yaev 
[19751 |LoveTace et al.|1999l ). 

Of course, the analytical formula (Eq. 1 ) describing the 
shape of the spot can be applied only in the cases where the 
magnetic field of the star has a strong dipolar component. In 
the case of a more complex field, the higher order component 
of the field may govern the funnel flow near the star, and 
the shapes of the spots will be determined by higher order 
multipoles of the field (e.g., Long et al.|2011 [[Romanova et al.| 
|2011| ) . This case should be investigated separately. 
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APPENDIX A: REFERENCE VALUES 

Our simulations are done using dimensionless variables (de- 
noted here by a tilde), obtained by dividing the dimensional 
variables by reference values. The reference values are deter- 
mined as follows: We choose fiducial values for stellar mass 
M.*, radius R* and equatorial surface magnetic field B*; the 
reference values of all other dynamical quantities are then 
obtained from these fiducial values. The unit of distance is 
chosen to be R = R*/0.35. The reference velocity is the Ke- 
plerian velocity at R , v = (GM,/R ) 1/2 , and uo = v /R is 
the reference angular velocity. The reference time is the Kep- 
lerian rotation period at R , Po = 27tR /v . 

The reference magnetic moment \i = where p is 
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CTTSs 


White dwarfs 


Neutron stars 




0.8 


1 


1.4 


R* 


2R 


5000 km 


10 km 


B* (G) 


10 3 


10 6 


10 9 


M- 


1 


1 


1 


Ro (cm) 


4 x 10 11 


1.4 x 10 9 


2.9 x 10 6 


Po 


1.8 days 


29 s 


2.2 ms 


B (G) 


43 


4.3 x 10 4 


4.3 x 10 7 


M (Moyr" 1 ) 


2.8 x 1CT 7 


1.9 x 10" 7 


2.9 x 10" 8 


F (erg s^cm -2 ) 


3 x 10 10 


5.6 x 10 17 


1.5 x 10 25 



Table Al. Sample reference values of the dynamical quantities used 
in our simulations in the case where dimensionless magnetic moment 
jl=l. Note that parameters po, Bo, Mq and Fq depend on jl (see 
Sec.[A). 



the dimensionless magnetic moment of the star. The reference 
magnetic field is B = M^/^o- The reference density is p = 
Bq/vq. The reference accretion rate is M = PoVoRo- Using the 
relationships for p , B and v , we obtain M = |a 2 /(p 2 v Ro). 

The dimensional accretion rate is M = MM , where M is 
the dimensionless accretion rate onto the surface of the star 
which is obtained from the simulations. Substituting the ex- 
pression for M and using \i = B*R 3 5 we get 



M : 



M\ B 2 R^ /2 (0.35) 7/2 



(GMJV2 



(Al) 



Thus, if the stellar parameters M.*, R* and B* are fixed, then 

the accretion rate is determined by M/p 2 . Therefore, vari- 
ation of parameter ft can be interpreted as variation of the 
magnetospheric radius of the star, or it can be equally validly 
interpreted as a parameter that is responsible for variation of 
the accretion rate. 

The reference energy flux is defined as F = PoVq. Us- 
ing this formula, we get a similar conversion formula for the 
energy flux: 




Figure B3. Dimensionless function F c = F c /jl 2 (which determines 
the energy flux in the centers of the spots) shown for all simulation 
runs. Symbols and lines are the same as is in Fig.|Bl| 



r m /R* are known, then Table Bl provides the positions and 
shapes of the spots. 

Figures [Bl] [B2] and [B3]s how parameters of the spots for 
all simulation runs. Fig. |B1| s hows the polar and azimuthal 
positions of the spots. Fi g. |B2| shows the polar and azimuthal 
widths of the spots. Fig. |B3| shows dimensionless energy flux 
in the center of the spots. A set of runs with the same 8 is 
marked by the same color. A set of runs with the same P* is 
marked by the same symbols. The lines show interpolations 
between points corresponding to sets of runs with the same 8 
and P*. These lines helped us derive important dependencies 
discussed in the paper (see Sec. [4}. 



,1/2 



(A2) 



We use a dimensionless parameter F c =Y c /\± 2 to charac- 
terize the flux in the centers of the spots. 



APPENDIX B: PARAMETERS OF SPOTS 

Parameters of the spots for all simulation runs are given in 
Table [Bl] Simulations were done for five misalignment angles 
8 = 2°, 5°, 10°, 20°, 30° and three (dimensionless) periods of 
the star P* = 1.8,2.4,2.8, which correspond to normalized 
corotation radii of r c /R* = 4.3,5.1,5.7. In each set of runs 
with the same parameters 8 and P*, we varied the dimen- 
sionless magnetic moment of the star jl and observed that 
the magnetospheric radius r m increases with p. The magne- 
tospheric radius is an important parameter because the fun- 
nel flow starts approximately at r m . In this paper, we use the 
dimensionless ratio r m /R*, so that results can be applied to 
different types of stars. 

Table |B1 1 shows the polar and azimuthal positions of the 
spots (do^o), their polar and azimuthal widths (9i, and 
the ratio F C /]I 2 , which characterizes the dimensionless energy 
flux in the centers of the spots. If parameters 8, r c /R* and 
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Table Bl. Parameters of the spots for different dipole misalignment angles 6, three periods of the star, P = 1.8, 2, 2.4, and different 
magnetospheric radii r m /R*. Polar and azimuthal angles (0 O and ty , respectively) determine the polar and azimuthal positions 
of the spots. Polar and azimuthal angles (0i and 4>i, respectively) determine the size of the spots. The dimensionless function 
F c = F C /|I 2 determines the energy flux in the centers of the spots. 
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Figure Bl. Left panel: The polar coordinates of the spots position O in all simulation runs as a function of the magnetospheric radius normalized 
to the stellar radius, r m /R*, for different misalignment angles: = 2° (red color), = 5° (green color), = 10° (blue color), = 20° (pink 
color), and = 30° (black color), and for different periods of the star: P* = 1.8 (thick solid lines, squares), P* = 2.4 (dashed lines, triangles), 
and P* = 2.8 (thin solid lines, circles). The symbols show results of individual simulation runs, while the lines show interpolation between runs 
performed for the same parameters and P* . 
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